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Abstract 

The Wheeler-DeWitt equation is obtained for Kasner-hke cosmologies. Some solu- 
tions to this equation are presented for empty space, space filled with a cosmological 
constant and in the presence of a scalar field. We also briefly discuss a non-commutative 
extension of these results. 
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1 Introduction 



The idea of quantizing gravity has been for long time being pursued and several approaches 
for this task have been developed. The canonical theory was one of the earlier attempts 
to follow this program and is based on a hamiltonian formulation of the field equations of 
general relativity |2l E] • One of the advantages of using a hamiltonian formulation is the 
fact that the notion of an absolute time does not exist, the time scale being arbitrary 
however the equations of motion do not suffer of ambiguities and they remain invariant under 
re-scaling. 

However several difficulties were present in this approach: it was not possible to associate 
conjugate momenta to some of the variables used in this formulation and therefore primary 
and secondary constraints had to be considered. The final outcome of this work were the 
dynamical constraints pn"^ = 0, H'^ = on a state vector The usual interpretation is 
that all the important information of the gravitational field that is coordinate-independent 
will be encoded on the hamiltonian constraint H'^ = 0. 

As in quantum mechanics, the wave function \E' is intended to determine completely 
the properties of physical systems of interest as quantum states evolving accordingly to the 
hamiltonian H. Due to the technical difficulty in treating the problem in all its generality, 
one is then led to analyze cosmologies that are distinguished because of the description of the 
universe that they might provide. In this context the Friedmann-Robertson- Walker universe 
has been studied in the presence of scalar and matter fields ^ . 

On the other hand it has been argued that near a cosmological time singularity the 
behaviour of a space-time is that of an ordinary manifold in empty space [H]. The study 
of homogeneous space-times with a group SU{2) of transformations acting on it has also 
provided some arguments in favor of this idea. One can also show that if any other fields 
are present they can be considered as small perturbations having as a background a vacuum 
solution. Because of the interest in knowing the evolution of the universe from its beginning 
as a singularity, the study of a Kasner-like cosmology using the quantization scheme provided 
by the canonical theory is then worthwhile. 

This paper is organized as follows: in section 2 we obtain the lagrangian and hamiltonian 
description for a Kasner-like cosmology. We use these results in section 3 to obtain the 
corresponding Wheeler-DeWitt equation and some solutions are presented in the case of 
empty space. Section 4 deals with the solution to the constraints p^"^ = 0, H"^ = in the 
presence of a cosmological constant. The case of a scalar field is then discussed in section 5. 
We end with some general remarks on the extension to the noncommutative setting of these 
results. 



2 Classical description 

We shall consider the following 4-dimensional Kasner-like metric 
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where represents the lapse function and the a^'s are functions depending only on t. Using 
the frame components 

5° = Ndt, 6' = e'^^dx' 

the Cartan formalism gives after a straightforward calculation the following expression for 
the curvature scalar 

R = -2N~^ [^{cui- oiiN-'^N + al) + ^ dtidij j . 

As in the usual way, the corresponding density Lagrangian from which we shall obtain the 
Wheeler-DeWitt equation is given by 



WttG 



N~'N + a'^) + J2'^iC^j) (2-1) 

i<j J 



where we have set SttG = 1. As discussed in [7], if one imposes A^ = 1 and afterwards one 
finds a solution to the variational equations 6L/Sai = 0, it is not sure that the action will 
be stationary at it. The appropriate way to proceed is to consider A^ as a free parameter 
and set A^ = 1 only after writing the dynamical constraints Pn"^ = and H'^ = 0. 

The second derivatives appearing in the previous expression can be ehminated by adding 
the total divergence 



div 



d_ 

dt 
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leading to the simple result 



OLi. 



(2.2) 



The associated momenta 



can readily be calculated, they are given by 



la N 



Using these results the Hamiltonian density can be written as 



N 



H = ^ aiTXa, - Lt = -— e ^ 



i<j 



Let us consider the constraint pn = dLx/dN = 0, which reads 

djdj = 0. 



i<j 



(2.3) 



(2.4) 
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As one can easily verify, the Kasner solution 

1 



a,; 



1 + a + 



1 + a, a(l + a), —a) Int + q = Pilnt + Ci 



(2.5) 



satisfies this constraint and thus the action we have just considered will be extremal at these 
values of the variables a^. 

The inclusion of a cosmological term is done without difficulty, following [7] one obtains 



H = N \ Ae^^ - - 



i<j 



and after integration over a space-like hypersurface t = const., one finds 

1 



H = NV \ Ae'^^ 
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2 J^TTiTT, 



i<j 



(2.6) 



where V is the associated 3-dimensional volume. Since this volume factor does not change 
at all the characteristics of the model, we shall disregard it in the next sections. We also 
have the following generalization of ()2.4j) 



aiCtj = A. 



(2.7) 



This condition can be obtained by considering again the constraint = dL/dN = 0, where 
now L is the lagrangian formed from ()2.2|) and L\ = —ANe'^. We shall consider some 
solutions to this equation in the following section. 



3 Quantization: Wheeler-DeWitt equation 

The standard procedure when quantizing a system described by a hamiltonian H consists 
in replacing the momenta vr by differential operators —iKV. This is a well known procedure 
from quantum mechanics, however, before doing that some words are in order. One of 
the ambiguities, if not the most important when translating quantities into their quantum 
counterparts is the fact that they do not longer commute. The Weyl prescription can be 
useful in this context, another possibility is to look for a kind of 'natural' ordering associated 
to the problem which is under study but a satisfactory answer to this does not exist. 

We shall chose a simple factor ordering and substitute each of the momenta appearing 
in (j2.6|) concordantly. This leads to the expression 



H 



i<j 



d 



(3.1) 



The simplest case one can consider for the above equation is that of empty space, i.e. 
A = 0. A solution to H'^ = is then being linear on the a^'s. This is a rather trivial 
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solution that is also unbounded on the space defined by these variables. It is then interesting 
to look for other solutions that can be regular. 



Consider the Ansatz 

V]/ = g«P"l+«'^"2+B(03) 

As a straightforward calculation shows, from the equation 



i<j 



^ = (3.2) 



it follows that the function -8(03) should satisfy the differential equation 

B + B^-2i{p + a)B-{p-af = 0. (3.3) 
Therefore, after a suitable translation on B one has 



B = i[p + 0") — u , u = —Aepa > (J. 



In order to further proceed we shall set Cq = ce'^, Ci = ce~^ with (p being interpreted as a 
constant phase. With this choice the fraction in the above expression becomes a hyperbolic 
function and a further integration gives then the final solution to (j3.3j) as 

5(0:3) = i{p + a)as - In | cosh.{pas + (f))\ + Bq. (3.4) 

Inserting this solution in the expression for ^ one has 

m = ^geip«i+i^a2+i(p+^)a3 / cosh(/ia3 + cj)) . (3.5) 

Therefore the probability density is given by 

I^P = |^olV(cosh(/ia3 + (t))f. (3.6) 

Some words are in order: first, were p negative then in (j3.4j) instead of cosh one would have 
cos and therefore the density probability l^^p would be divergent at the points 03 satisfying 
/ia3 + = (2n + l)7r/2, with n an integer. The fact that ()3.6p is bounded and vanishing when 
ai — i> ±00 points out to the possibility that regularized solutions can exist for more general 
situations other than empty space. Another important point concerns the dependence on 
03. As can be seen from fKllj) . the Hamiltonian is invariant under a cyclic permutation of 
the indexes 1, 2 and 3. It follows then that expressions similar to ()3.(jj) with ai or 02 instead 
of 03 exist. 

It is a straightforward calculation to find the current J = y(\E'(9\E'* — \E'*(9\E') associated 
to (j3.5j) . One has 

J = /i|^P(p, (T,p + a) 

which is also bounded. 
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4 Cosmological constant 



Before solving the Wheeler-DeWitt equation ()3.ip in more generality, we shall consider the 
equation of motion for the variables a^, namely ()2.7|1 . In order to solve it we take the 
following Ansatz 

Fit) 

ai{t)=pM-=;f-\ + qiK{t) 

with the Pi's and g^'s being for the moment arbitrary constants. Here to is an arbitrary time 
to be specified later on and that plays the role of a scaling factor. Using this expression one 
arrives to the equation 



A 



(4.1) 



where the quantities 



P' 



have been defined. 



i<j 



i<j 



i<j 



We know that in the case of a vanishing cosmological constant the constants pi can be 
chosen to be the Kasner exponents once one sets F{t) = Int and K{t) = 1 (see (j2.5p ). In 
the following we would like to consider the departures from this "classical" solution for small 
values of the cosmological constant. For this we shall solve K{t) in terms of F{t) using ()4.1|) . 

There are two different cases to be considered depending on the value of q^: 

a) q^ = 

There is not quadratic term on K in ()4.1|) and one can write directly 



Fit, 



^dt. 



to 



The solution is unique and one then has 

= (F(t)/F(to))^^"'''^exp 



Kit, 



qiA 



-^dt 



to 



(4.2) 



In this case there are two different solutions 
Accordingly we have 



to 



A m'^ - Ap^q^ 



4g4 F2 



1/2 



dt. 



e"^± = (F(t)/F(to))"^exp 



K^ito)±q^ 



to 



A ^ - Ap^q^ F^ 



q^ 



4g4 p2 



1/2 



dt 



(4.3) 
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Let us now specialize these solutions to the Kasner case, namely = and F{t) = t. 
We also choose to = 1 so that -F(to) = 1- Case a) gives then 



2 m? 



By setting K{1) = A/2m^ one has for A ^ 1 



2 

up to linear terms. It follows that to lowest order the volume element is then given by 

1 rA 



dVia) = V^d^x = Nt{l + -—t 

2 m'^ 



(4.4) 



In the above expression we have used the fact that 'Y^^Pi = 1. 
For case b) one has 



KM) 



2 

m 

V 


Int ± 


2 

m 


2 

m 




2 

m 


lnt± 



1 + 



4g2A 



1/2 



dt 



1 + (1 + 



2^/2 + ln(l + 2' 



where = Aq^ At^ / . For small A we have 

1 A 



2 m'' 2 



where we have chosen 



K+(l) = = -^[1 + In \qA'/ym\ - 2'/^ + ln(l + 2^^% 

2q 



Using the above results one obtains 

2 



2 



and therefore the corresponding volume elements are 



1 rA 

dV(b)+ = Nt(i + -—e 



2 m? " 



1 rm^ 1 rA r, 



2m^ 



(4.5) 



We see from ()4.4j) and ()4.5|1 that to lowest order on A, the volume elements dV(a) and 
dV(^b)+ are equal while dV^by has a supplementary dependence on t (the values of K{1), K±{1) 
have been chosen as to eliminate all the irrelevant constant scaling factors from these ex- 
pressions). We can now define the ratio rji = dV(^i)/dVK, with dVx = Nt the volume element 
associated to the Kasner cosmology, in order to compare the different geometries. We have 

1 rA r, -rs^ , 1 rA 9, , , 
= = 1 + -—t\ r]b-=t ^{1- -—t"). 4.6 

2 m'^ 2 
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It is clear that different behaviors in the hmit t — are present: in the first case one can 
effectively replace dV(^a) or dV(^b)+ by dYx- The second case is more interesting since due to 
the presence of the exponent z = rm? /q^ the ratio can either vanish {z < 0) or diverge 
(z>0). 

We shall now turn our attention to After a translation ctj — > + q with /i^/4 = 

Ae^^'^'/S, it becomes 



^{aA = 0. 



(4.7) 



A solution can be obtained by considering \1/ = \1/(C)- Since Vj = V,^ one arrives to the 
differential equation 

The equation can be solved by setting x = e^. Using this variable one obtains a Bessel 
differential equation of order 0. The solution is then simply given by 



^ = AJo{x) + BYo{x) 



(4i 



with Jq{x),Yq{x) the Bessel and Neumann functions of order respectively. We have the 
following asymptotic developments 



\l/(x) ~ —;={A cos{x — 7t/A) + B sin(x — 7r/4)), x — > oo 



and 

^(a;) ~ A + 51n(x/2), x 0. 

As in the case of empty space, a bounded solution for all values of C {x & [0, oo)) is possible 
to exist: it suffices to take A ^ 0, B = 0. 



5 Inclusion of a scalar field 



When massless scalar and vector fields are present, it is known that the behavior of the 
Kasner exponents in the classical solution is modified [H]. Through evolution in time, the 
classical oscillatory regime, known as Kasner epochs [Hj and characterized by having one of 
the Pi's to be negative is modified and in fact, a final Kasner epoch is reached where all 
the Pi's take positive values. The reason for this is that in the early stages of the cosmic 
evolution, a massless scalar field will give a contribution not negligible to the total stress- 
energy tensor, roughly of the same order as the gravitational interaction. 

The kinetic energy associated to a scalar field $ depending only on t is given by 
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and therefore its contribution to the quantized hamiltonian ()3.1|) reads 



Hg, = —III 



2 * 2 

The full Wheeler-DeWitt equation to be considered is then 



2 



V| 



+ V{^)e^^ ^ ^ = 0. 



-i<3 ■■ 

We rewrite the first two terms in the brackets with the help of a metric g^^ as 



(5.1) 



i<j 




Through the usual trick of diagonalization one has 



Xl 



+ 1/(<l>)e2^|^ = 



(5.2) 



where now 



C = ^[(1 + + (1 + ^)^2 + (1 - 2x)x3], 

V3 



X 



We have also performed a similar translation as in ()4.7|) in order to eliminate a factor (for 
details see the Appendix). 

The above transformation is suitable to find some solutions to (j5.2p both for V{^) = 
and V{^) = A. We shall not dwell on this, instead we shall concentrate on the solutions to 
the classical equations of motion in the presence of the scalar field In order to do so we 
shall stick to the variables since they are better adapted to make the connection with the 
Kasner metric. We have accordingly 



i<j 



We remark that the problem contains now four unknown functions to be determined, however 
we only have at our disposal two equations. Instead of trying to find the most general solution 
we shall take a similar Ansatz as that of section E] 

ai{t) =pi\nt + qiK{t) (5.3) 

with Pi being a Kasner exponent. We shall also set V($) = A. By doing this one arrives to 



mHk + qH^k^ = ^e-^'^'' + At^, 



Making the change of variable s = e one finally obtains 



m 
2r 



47-2 ' 



'-^tss + = —s + Kt^s^. 



(5.4) 
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In this last equation we have also supposed r 7^ (otherwise with some little modifications 
we can use directly the results of the precedent section to write down a solution). 

Unfortunately we have been not able to find the most general solution to ()5.4|) . However, 
for the special case = one obtains 

s{t) = cie'-''*'/"' - I^e'"''"/"'" Ei(l, rAf/m^). (5.5) 

The function Ei(n, x) can be written as an infinite series on x and is ill-defined at x = 0, it 
has a logarithmic divergence. Due to the fact that $ ~ t^^s^^^"^, it follows also that for this 
solution the scalar field $ can not be defined for values of t smaller than a cut-off > 0. 



6 Conclusions 

We have found the Wheeler-DeWitt equation for 4-dimensional Kasner-like space-times and 
some cases of interest have been considered. For empty space the main result is that well- 
defined solutions exist in the space of the parameters of the model and the equations of 
motion admit the well known Kasner metric as solution. When a potential corresponding 
to a cosmological constant is introduced, a simple quantum state \l/ can be found which has 
finite density probability. The solution to the classical equations of motion for this case show 
that different behaviors for the geometry of the space are present in the limit t — *• 0. With a 
proper choice of the parameters involved the volume element of each solution can be made 
bigger than that of the Kasner metric. In particular when < z < 1, the ratio r]i,_ diverges. 

The presence of a scalar field presents some additional complications. However, to write 
down a explicit solution should not present any special difficulties using equation ()5.2|1 . On 
the other hand, the classical equations of motion allow for a solution for the scalar field in a 
simplified case. However the solution thus found is defined only in a certain region {t > to) 
and shows the kind of problems one may encounter in the most general case. 

A generalization of the above results to dimensions greater than 4 should be possible in a 
straightforward way. We would like to end with some comments on a non-commutative exten- 
sion of these results. The interest in introducing a non- commutative structure for 'classical' 
theories in some respect comes from the natural generalization it provides. From the physical 
point of view this has been important in giving new insight to some old problems pilj illj . 
Here we are mainly interested in the link with gravity that can be stablished ^2] and the 
natural cut-off, essentially of the order of Planck's length, that it might provide to regularize 
singularities either space or time-like [T^ I14j . 

Examples of the role non-commutativity can play in cosmology are to be found in 
EI- The most natural way to introduce non-commutativity is through the use of the 
Moyal product, which in its differential form is given by |18j 



= exp 



2 " 



fix)gix')\,-. 
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The non-commutative parameter 9 is defined through the commutation relations 
iOij. With this one can then write 



i<j 



where Vt is an arbitrary potentiaL The non-commutativity of the a^'s can be translated into 
a redefinition of these variables by making the replacement — >■ + ^OijTi^ . This change 
will affect only the term Vre^^. An account of the consequences of this will be presented 
elsewhere. 



Appendix 

From the metric g"^^ introduced in section El we have 

/O 1 1\ 

9AB=\l 1 . 

V 1 V 

The eigenvalues of qab are —1 (multiplicity 2) and +2. The matrix diagonalizing g^s is 

^ 1 2x —X —x\ 
D-^ = — I -X 2x -X , x = J2/3 
^[1 1 ij 

and hence qab = {D~^gD)AB = diag(— 1, —1,2). It follows that 

Va = -QabV^ = ~gAB{D~')f,g'^^'yM = {D-')iVB. 

Explicitly we have 

v., = ^[2Vi-V2-V3] 

V,. = ^[-Vi + 2V2- V3] (A.l) 
Vx3 = ^[Vi+ V2 + V3] 

and 

xi = -^[ai + xa3] X2 = -^[ai + xa3\ X3 = ^[-ai - a2 + xa^]. (A. 2) 
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